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Abstract. We fit WMAP5 and related data by allowing for a CDM-DE coupling and non-zero neutrino masses, simulta- 
neously. We find a significant correlation between these parameters, so that simultaneous higher coupling and v-masses are 
allowed. Furthermore, models with a significant coupling and v-mass are statistically favoured in respect to a cosmology 
with no coupling and negligible neutrino mass (our best fits are: C ~ l/2mp, ~ 0.12eV per flavor). We use a standard 
Monte Carlo Markov Chain approach, by assuming DE to be a scalar field self-interacting through Ratra-Peebles or SUGRA 
potentials. 

Keywords: dynamical dark energy, massive neutrinos 
PACS: 14.60.Pq,95.36.+x,95.35.+d 

WHY DE AND CDM SHOULD BE COUPLED 

One of the main puzzles of cosmology is why a model as ACDM, implying so many conceptual problems, apparently 
fits all linear data OlElsIl in such unrivalled fashion. 

It is then important that the fine tuning paradox of ACDM is eased in cosmologies where DE is a self-interacting 
scalar field (j) (dDE cosmologies), with no likelihood downgrade |4]. The coincidence paradox is also eased in coupled 
DE (cDE) cosmologies, i.e., if an energy flow from CDM to dDE occurs |5]. CDM-DE coupling, however, cuts the 
model likelihood, when we approach a coupling intensity significantly attenuating the coincidence paradox. 

The physical cosmology could however include a further ingredient, able to compensate coupling distorsions. Here 
we show that a possible option is neutrino mass. In fact, when we assume CDM-DE coupling or a significant v mass, 
we cause opposite spectral shifts f^. 

It is then natural to try to compensate them; if we do so, however, the residual tiny distorsions tend to favor coupling 
and V mass, in respect to dDE or ACDM. 

We illustrate this fact by considering the self-interaction potentials 

y((^))=A«+V(^)« (1) 

or 

Vi^) = (A«+V</)«)exp(47r(^)Vm2), (2) 

(RP fl^ and SUGRA |8], respectively; nip : the Planck mass), admitting tracker solutions. Uncoupled RP (SUGRA) 
yields a slowly (fastly) varying w{a) state parameter. Coupling is however an essential feature and modifies these 
behaviors, mostly lowering w{a) at low z, and boosting it up to +1, for z >^ 10. 

For any choice of A and a these cosmologies have a precise DE density parameter Here we take A and as 
free parameters in flat cosmologies; the related a value then follows. 

In these scenarios, DE energy density and pressure read 

P=Pa- + V(0), p^Pk^V{(j)), (3) 



with 



Pk - f/2a^ ; 



(4) 
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FIGURE 1. Left panel: Transfer functions (multiplied by Z:''^, for grapliic aims) in example cosmologies with/without coupling 
and witli/witliout 2 massive v's (00/CM models). Right panel: Binned anisotropy spectral data, normalized to the best-fit ACDM 
model (upper frame) or to the best fit SUGRA cDE model including v-masses (lower frame). The distorsions arising when coupling 
or v-mass are separately considered are also shown. 



dots indicating differentiation in respect to T, the background metrics being 

ds^ ^a'^{'z)[d'Z^ -dX^] (5) 

with 

dX^ ^dr^ + r^{dd^ + cos^dd<^^) . (6) 

Until ^ V , therefore, w{a) — p/p approaches +1 and DE density would rapidy dilute (p °^ a^^), unless a feeding 
from CDM occurs. The gradual increase of (j) leads it to approach nip, so that aV ^ Pk regime is attained. The state 
parameter approaches then -1 and DE induces cosmic acceleration. 

DE cannot couple to baryons, because of the equivalence principle (see, e.g. J^). Constraints to CDM-DE inter- 
actions, however, can only derive from cosmological data. The simplest coupling is a linear one, formally obtainable 
through a conformal transformation of Brans-Dicke theory (see, e.g., 1 10]). The coupling intensity must however be 
adequate to transfer from CDM to DE the energy needed to beat its spontaneous dilution a^^. This prescribes a 
coupling scale C^' ^ nip, while the energy drain from CDM lets its density decline (slightly) more rapidly than a^^. 
The whole scenario between recombination and the birth of non linearities is then modified, and it comes as no sur- 
prise that we expect significant changes both in the matter density fluctuation spectra P{k) and in the observed angular 
anisotropy spectrum Q. 

The prescription that the 2-component dark sector interacts with baryons or radiation only through gravitaion reads 

T^'^K'^li + 7^' ''' v(Ji = (here Tjly ''^^ are the stress-energy tensors for CDM and DE, let their traces read r''^''''''). Only 
if we make the further assumption that the 2 components are separate, it shall be C = in the relations 

T^''%i^^ = +CT^'^^^ , - -crM^.v , (7) 

obtained by passing from the Jordan frame, where a Brans-Dicke frame cosmology holds, to the Einstein frame ifsi fToll . 
These equations also show why DE cannot couple to any component with vanishing stress-energy tensor trace, e.g. to 
radiation. 

Besides of C, we shall also use the dimensionless coupling parameter p — (3/16;r)'/^mpC . If it is C = l/nip, 
therefore, )3 = 0.2443; values )3 ^^{0.1) mean then C ~ 1 /Inip . 
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FIGURE 2. Marginalized (solid line) and average (dotted line) likelihood of cosmological parameters in SUGRA models. Notice 
that the /3 signal appears when high- and low-z data are put together, and is strengthened by SNIa data. As a matter of fact, coupling 
allows to lower the "tension" between flc and fluctuation amplitude detected from CMB and deep sample data. For My we also 
show the corresponding likelihood distributions obtained in the case of a standard ACDM+My model (gray lines). 



WMAP5 




0.1 0.2 0.3 -10 10 0.5 1 1.5 2 2.5 



WMAP5+2dF 




0.1 0.2 0.3 -10 10 0.5 1 1.5 2 2.5 



all data 




0.1 0.2 0.3 -10 10 0.5 1 1.5 2 2.5 

3 log^p(A/GeV) MyeV 



FIGURE 3. As previous Figure, in RP models. 

In this context let us then remind the interesting option of considering DE coupled to v's ifTlll ; the trace T^^^ becomes 
significant only when their kinetic energy is redshifted to values ^ niy. 

Another option considered in the literature is that the r.h.s.\ of eqs. (|7]l are replaced by iCT '''''' 0,v, where C sign 
is opposite to C lfT2ll . This option causes a bootstrap effect, as the energy transfer to DE is boosted as soon as its 
dilution is attenuated. In analogy with the coupling to v's, this approach therefore associates the peculiarity of our 
epoch (the coincidence) with another accepted peculiarity. Furthermore, at variance from our approach, this option 
does not follow from any conformal transformation of Brans-Dicke formulation. 

It should be however noticed that the correlation between My (the trace of the 3-dimensional neutrino mass matrix) 
and C seems no longer to persist within the frame of this approach if new datasets are considered 1, 1 3.1 . 




FIGURE 4. Two parameter contours for the SUGRA model. Solid lines are 1- and 2-(T limits for marginalized likelihood. Colors 
refer to average likelihood, and the 50% likelihood contour from the average likelihood is indicated by the dotted line. 




FIGURE 5. As previous Figure, for a RP potential. 



PROCEDURE AND RESULTS 

The results shown in this paper are based on our generalization of the public program CAMB 1 14], enabling it to study 
cDE models. Likelihood distributions are then worked out by using CosmoMC 1 15]. 

In Figure [T] we then illustrate the compensating effects between coupling and v mass, by considering various 
cosmologies. A 00-model is then a cosmology without coupling and v mass. In turn, the 01- and 10-models are 
cosmologies with coupling (and no v mass) and v mass (and 0-coupling). In the 1 1-model, finally, both v mass and 
coupling are included. 

The l.h.s. plot shows an example of transfer functions for 00-, 01-, 10- and 11-models. The sum of v masses 
(My) is tuned so that the 00- and 1 1- r(fc)'s nearly coincide. At the r.h.s., we then see a similar comparison with C/ 
anisotropy spectral data. More details are provided in the caption. 

More details on the fitting procedure are provided elsewhere [17]. The fitting procedure returns parameter values 
mostly in the same range as for dDE or ACDM cosmologies. The significant parameters in our approach are however 
the energy scale in potentials (A), the coupling intensity ()3), and the sum of v masses (My). In Figures |2] and [3] we 
provide one-dimensional hkeUhood distributions on these parameters for SUGRA and RP cosmologies. 
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FIGURE 6. Evolution of density parameters in a SUGRA model with coupling and v mass. Colors refer to different components, 
as specified in tlie frame; lines to the different models: My =0, j5 — (continuous line); My = 0.5 eV, j8 = 0.085 (dotted); 
Afv = 1.1 eV, j8 = 0.17 (short dashed); My = 1.2 eV, j3 = 0.22 (long dashed). 
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FIGURE 7. Detection of CDM-DE coupling, following an hypothetical determination of an electron neutrino mass my — 0.3 eV, 
by the experiment KATRIN, yielding My ~ 0.9 eV. Here we show the likelihood distribution, with such a prior on My. Colors and 
lines convey the same indications as in Figures 4 and 5. 

A basic information is however the correlation between likelihood distributions. These are shown in Figures|4]and 
|5]again for SUGRA and RP cosmologies, respectively. 

These Figures, as well as one-dimensional plots, clearly exhibit maxima, both for average and marginalized 
likelihood, for significantly non-zero coupling and v-masses. Although their statististical significance is not enough 
to indicate any "detection" level, the indication is impressive. Furthermore, a stronger signal, with the present 
observational sensitivity, would be impossible. 

Let us however outline that this work aimed at finding how far one could go from ACDM, adding non-zero coupling 
and V-mass, without facing a likelihood degrade. It came then as an unexpected bonus that likelihood does not peak 
on the 0-0 option. 



CONCLUSIONS 



The allowed j3 values open the possibility of a critically modified DE behavior. Figure |6] shows the scale dependence 
of the cosmic components for various jS-My pairs. 

For P values comprised between 0.1 and 0.22, i.e. well within 2-(7's from the best-fit model, just as the 0-0 option, 
we have a long plateau in the energy density of DE, going from z ~ 10 to above z ~ 1000. The devline of p^/g at greater 
z is then mostly due to a parallel decline of the density Pc of CDM: when radiation becomes dominant, both DCM and 
DE are negligible. 

The ratio Pc/pde in the plateau is i^(lOO). At lower z it becomes gradually smaller because of the gradual 
contribution of the potential term to p^^. It should be however outlined that this behavior is not ad-hoc, but is incribed 
in the tracker solutions for the potentials selected. 

Let us finally outline that the My values allowed by j3 ^ 0.1-0.2 approach the V-mass detection area in the 
forthcoming experiment KATRIN llT6ll . based on tritium decay. 

Accordingly, should particle data lead to an external prior on My, the strong degeneracy between the coupling 
parameter )3 and the neutrino mass My is broken, and new insight into the nature of DE is gained 1 17]. In Figure|2]we 
show how a neutrino mass determination symultaneously implies an almost model independent CDM-DE coupling 
detection. 

This would be a revival of mixed DM models ifisll . in the form of Mildly Mixed Coupled (MMC) cosmologies. 
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